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Abstract 

Relationships between diverse thermodynamic quantities appropriate to weakly compressible 
turbulence are derived. It is shown that for turbulence of a finite turbulent Mach number 
there is a finite effect of compressibility. A methodology for generating initial conditions 
for the fluctuating pressure, density and dilatational velocity is given which is consistent 
with finite Mach number effects. Use of these initial conditions gives rise to a smooth 
development of the flow, in contrast to cases in which these fields are specified arbitrarily or 
set to zero. Comparisons of the effect of different types of initial conditions are made using 
direct numerical simulation of decaying isotropic turbulence. 
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1. Introduction 

Compressible shear layers are encountered in many practical applications ranging from su- 
personic ejectors for mixing and entrainment in gas turbines to scramjet combustion. Such 
flows, important in mixing enhancement and noise reduction, are modified using different 
nozzle shapes or strategies such as multiple jets, coaxial jets, and countercurrent mixing 
lavers. Compressible mixing layers, for a certain class of initial conditions, have a growth 
rate substantially lower than their incompressible counterparts. Compressibility effects also 
contribute to substantial reductions in turbulence intensity levels and the turbulent shear 
stress anisotropy while increasing levels of anisotropy of the normal stresses. Additional 
details regarding the effects of compressibility can be found in Bradshaw [1], Lele [2], Spina 
et ai [3], and Gutmark et al. [4]. 

Many of these flows are part of a class that have a small bulk dilatation and in which 
turbulence gradients are more important than mean pressure gradients. Such flows typically 
have a small turbulent Mach number squared, M ( 2 <C 1, where M t = u/c, u is a turbulent 
velocity scale and c is the mean speed of sound. This reflects about the weakly compressible 
nature of the turbulence — the mean flow itself may be highly supersonic. An analytical 
development of the covariance with the fluctuating dilatation suitable for such a class of 
flows has been given in Ristorcelli [5]. This note explores the implications of that analysis 
in the context of the initial conditions for DNS or turbulence model calculations. It is 
shown that as a consequence of a finite turbulent Mach number there are finite non-zero 
density, temperature and dilatational fluctuations. Moreover these fluctuations are, to lowest 
order, specified by the incompressible field. These results are relevant to initial conditions 
on thermodynamic quantities employed in DNS of compressible turbulence that are started 
using incompressible fluctuating velocity fields. A potential consequence of the present results 
is a reduction in the amount of computational effort spent adjusting to transients associated 
with the relaxation from arbitrary initial conditions. Also, for flows such as decaying isotropic 
turbulence which are highly dependent on initial conditions, the methodology presented here 
provides a more natural type of initial conditions than others which have been previously 
considered. 

2. Analysis 

Ristorcelli [5] has conducted a small Mach number expansion of the compressible Navier- 
Stokes equations. Some of those results, as they are relevant to the present subject, are 
now briefly sketched. The problem of compressible turbulence can be viewed as a singular 
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perturbation. The inner problem is related to the acoustic source problem with time and 
length scales based on the turbulence time and length scales, lju and t. Here u and £ are 
a fluctuating velocity scale and a turbulence integral length scale. The outer problem is 
the acoustic propagation problem, which in the low M 2 limit has time and length scales 
A/c = C/u and A = £/M t . The inner expansion is relevant to the compressible turbulence 
modeling problem for compact flows [5]. 

The velocity field is decomposed according to u, = v t + e 2 u>; + ... thus d = u,,i = c 2 w t ,i +..., 
where the small parameter is e 2 = 7 M 2 where M t = u/c. The thermodynamic variables are 
decomposed according to a perturbation about a mean state, (P,p,T), thus p* = P( 1 +p), 
p* = p(\ + p) and T* = T{ 1 + 9). Then ansatzs of the form p = e 2 [pi + e 2 p 2 + •••] are used. 
In the context of turbulence models, the following definition of the fluctuating Mach number 
is used by M 2 = §fc/c 2 where k = 1/2 < UjUj >. To lowest order 
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In light of the homogeneous shear flow simulations of Blaisdell et al. [ 6 , 7] (the fluctuations 
were seen to follow a polytropic gas law with coefficient close to the adiabatic value), the 
adiabatic case is treated. The dimensional form of the linearized adiabatic equation of state 
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accounting for the last two equations in the zeroth-order expansion above. The zeroth-order 
problem is the incompressible problem. The fluctuating pressure is not an independent 
variable but is set by the vortical velocity fluctuations and which also produces the density 
and temperature fluctuations. On the inner scales the incompressible pressure fluctuations 
dominate the pressure field. This pressure is called by acousticians the “pseudo-pressure" 
in contradistinction to the propagating pressure associated with the sound field. The next 
order expansion for the continuity equation produces a more complex diagnostic relation for 
the fluctuating dilatation, 
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The relationships between the density, temperature and dilatation as well as the velocity 
field associated with the dilatation to be derived presently will be called “pseudo-sound ,? 
relationships. The dimensional equivalents of the above relationships, using the facts that 
pi = ffjp, 81 = O'/T. pi = p'/P are 
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where p' is the incompressible pressure fluctuation. In the context of a calculation two 
different nondimensionalizations are possible. In either case the pressure fluctuations will 
scale with pu 2 . For Case 1 the density and temperature fluctuations are scaled with the 
turbulence fluctuations, pu 2 /c 2 and u 2 /c p respectively, and one obtains p = p and 6 = p. 
For Case 2 the density and temperature fluctuations are scaled with the mean temperature, 
T, and density, p , and one obtains p = M 2 p and 6 = (7 — 1 )M 2 p. In either scaling 
the fluctuating dilatation is obtained from -d = M?(p, t +v k p, k ). In both cases it is clear 
that for finite non-zero turbulent Mach number there are finite non zero fluctuations of 
density, temperature, dilatation and pressure and there are specific relationships between 
these quantities. 

Obtaining an initial condition for the compressible portion of the velocity field, in,, which is 
consistent with the relation (7) is somewhat involved. The fluctuating pressure pi needed in 
(7) is obtained from the Poisson equation (3), and its time derivative can be found by taking 
the time derivative of (3) and then substituting for the time derivative of the velocity from 
(1). This give a Poisson equation for pi, t , 

(piit ) ijj — 2 [{v k Vi, k "Fpi,: )uj] ,jj (11) 

Solving the two Poisson equations (3) and (11), the dilatation is then found from (7). For 
homogeneous turbulence the dilatational velocity can be found from the dilatation by working 
in Fourier wave space. Since the dilatational velocity is irrotational the Fourier coefficients 
of the velocity are aligned with the wavevector, k , and they can be found from 

Wj = ( 12 ) 

The dilatational velocity tc; is then combined with the solenoidal velocity u, to obtain the 
full velocity field. 
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3. Simulations and Results 


In order to determine the effectiveness of the new initial condition method, direct numerical 
simulations of decaying isotropic turbulence are done with three different types of initial 
conditions. For all three types of initial conditions a solenoidal velocity field is specified with 
an initial three-dimensional spectrum given by 

E(k) = Ak 4 exp(—2k 2 /k 2 ) (13) 

and random phases for the Fourier coefficients. The constant A determines the amplitude of 
the velocity fluctuations and is set to obtain a specific turbulent Mach number, M t . The peak 
in the spectrum is given by the wavenumber k p which is set to 12 for all of the simulations 
presented here. The first method, ICl, specifies zero density and pressure fluctuations. The 
velocity field is purely solenoidal so that the dilatation is zero. The second method, IC2, has 
zero density fluctuations, but the fluctuating pressure field is found by solving the Poisson 
equation (3). The velocity field is again purely solenoidal. The third method, IC3, uses the 
methodology developed here. The fluctuating pressure field is found by solving the Poisson 
equation (3). The density fluctuations are related to the pressure through (4). The dilatation 
is nonzero and is found from (7). Lastly, the dilatation is used to find the dilatational velocity 
from (12). 

Details about the simulation method can be found in Blaisdell et al. [6]. Simulations were 
performed with initial turbulent Mach numbers M t = 0.231, 0.115 and 0.058. The initial 
turbulent Reynolds number Rej = pq 4 /ev = 200 where q 2 = pu t u x /p and e is the dissipation 
rate of turbulent kinetic energy per unit volume. 

For the present brief communication one figure, using data from the M t = 0.231 runs, has 
been chosen to illustrate the analysis. Figure 1 shows the development of the variance of the 
dilatation, Id. It is clear from the figure that there are sizable time boundary layers near the 
origin for the first two methods as the flow seeks to adjust to inconsistent initial conditions. 
This is not the case for the pseudo-sound initial conditions. Also interesting to note is the 
fact that the variance of the dilatation for ICl is one order of magnitude larger than IC3. In 
general, at lower M t this difference is more pronounced and at M t = 0.058 the difference is 
two orders of magnitude. This excess dilatation is associated with the initial condition and 
not the dilation due to the evolving compressible velocity field. These results are consistent 
with the physics of the near field evolution as developed in [5]. 
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An intuitive argument can be given to explain the behavior seen in Fig. 1. One might spec- 
ulate that there is a wave field generated by the inconsistency between the pressure, density, 
temperature, and dilatational fields whose decay rate is far slower than the turbulence decay. 
This gives rise to a background acoustic radiation that lasts the course of the simulation. 
For IC1 the fluctuating pressure field is set to zero. This can be viewed as a combination 
of the incompressible pressure satisfying (3) plus an acoustic pressure field which exactly 
cancels the incompressible pressure field. As time evolves the two pressure fields become 
decorrelated, giving rise to large acoustic pressure and dilatation fluctuations. For IC2 the 
pressure field is correct; however, there is no dilatation field. Again this can be viewed as 
a combination of the dilatation found from (7) plus an acoustic dilatation field that exactly 
cancels this. As the flow evolves away from the initial conditions, the two dilatation fields 
become decorrelated so that one is left with the correct dilatation plus acoustic fluctua- 
tions of dilatation and pressure. The evolution of the pressure variance (not shown here) 
corroborates this picture. 

4. Conclusions 

The major result of this note is that the initial conditions on the fluctuating thermody- 
namic variables of density, temperature, pressure, and dilatation should not be arbitrarily 
specified in the DNS of compressible turbulence. To lowest order these fluctuations, whose 
nature is connected by the underlying fluctuating vortical turbulence field, are generated by 
the pressure field associated with the divergence-free portion of the vortical motions. The 
methodology for generating initial conditions presented here allows the flow to develop more 
naturally. The significance of the initial conditions may depend on the type of flow consid- 
ered; however, one is led naturally to speculate about the nature and relevance of simulations 
that start from arbitrary initial conditions. 

A discussion of initial condition issues is given in Blaisdell et al. [7]. In homogeneous shear 
flow the effects of arbitrary initial conditions are seen to dissipate by St = 6, where S is the 
shear rate. In isotropic turbulence, however, the initial conditions do not fade. Blaisdell et 
al. [7] has pointed out that the lack of coupling between the wave and vortical fields in the 
isotropic flow is a probable source of these two different dependencies on initial conditions. 
(This should also be understood as an additional dependence on the initial conditions other 
than that already established for incompressible isotropic turbulence.) 

There are surely a large number of interesting studies in which one can study the relaxation 
from diverse arbitrary initial conditions. It is also possible to argue that they have practical 
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value. However, the strongest argument, given the current engineering problems and the lack 
of knowledge regarding the effects of compressibility, can be made for the initial conditions in 
which the density, temperature, pressure, and dilatational fields are related to the underlying 
local turbulence field that is, to a large degree, the source of the fluctuations. This seems 
better than an arbitrary guess at initial conditions that may or may not have physical 
relevance or meaning and may produce unrealizable results. 
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